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Systematic description of a spin one-half system endowed with magnetic moment or any other 
two-level system (qubit) interacting with the quantized electromagnetic field is developed. This 
description exploits a close analogy between a two-level system and the Dirac electron that comes 
to light when the two-level system is described within the formalism of second quantization in terms 
of fermionic creation and annihilation operators. The analogy enables one to introduce all the pow- 
erful tools of relativistic QED (albeit in a greatly simplified form). The Feynman diagrams and 
the propagators turn out to be very useful. In particular, the QED concept of the vacuum polar- 
ization finds its close counterpart in the photon scattering off a two level-system leading via the 
linear response theory to the general formulas for the atomic polarizability and the dynamic single 
spin susceptibility. To illustrate the usefulness of these methods, we calculate the polarizability and 
susceptibility up to the fourth order of perturbation theory. These ab initio calculations resolve 
some ambiguities concerning the sign prescription and the optical damping that arise in the phe- 
nomenological treatment. We also show that the methods used to study two- level systems (qubits) 
can be extended to many-level systems (qudits). As an example, we describe the interaction with 
the quantized electromagnetic field of an atom with four relevant states: one S state and three 
degenerate P states. 

PACS numbers: 12.20.Ds,42.50.Ct,32.80.-t,76.20.-|-q 



I. INTRODUCTION 
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(b) 



Two-level quantum systems, called qubits by Schu- 
macher |l| , play a fundamental role in quantum informa- 
tion theory. In this context they are usually treated as 
mathematical objects living in a two-dimensional Hilbert 
space. In reality, qubits always exist as material objects 
and we should not forget that they are endowed with con- 
crete physical properties. In this paper we shall deal with 
two-level systems that interact directly with the electro- 
magnetic field, such as spin one-half particles endowed 
with magnetic moment or two-level atoms. Thus, our 
results do not apply to qubits encoded in the polariza- 
tion states of photons. We shall restrict ourselves in this 
paper to isolated qubits interacting only with the quan- 
tized electromagnetic field. Therefore, the calculated de- 
cay rates will include only the spontaneous emission. 

A two-level system is the simplest model of a quantum 
system and yet in the presence of a coupling to the quan- 
tized electromagnetic field an exact solution has not been 
obtained. Even in the simplest case, when the electro- 
magnetic field is restricted to just one mode, the model 
has been exactly solved only in the rotating-wave ap- 
proximation by Jaynes and Cummings Among the 
approximate solutions, perturbation theory is still the 
most universal and effective tool, especially in the world 
of electromagnetic phenomena. 

In the present paper we develop a systematic and com- 
plete theory based on an observation that a two-level sys- 
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FIG. 1: Two Feynman diagrams representing the elementary 
processes and their interpretation in terms of the Dirac-sea 
picture. The pair creation (a) corresponds to the photon ab- 
sorption causing a transition (b) of the two-level system from 
the ground state to the excited state. The electron is moved 
from the negative energy state (creating a hole) to the pos- 
itive energy state. The pair annihilation (c) corresponds to 
the inverse process (d). The electron jumps back from the 
positive to the negative energy state emitting a photon. 



tem can be treated as a relativistic trapped electron. The 
translational degrees of freedom of such an electron are 
practically frozen. The only "degree of freedom" that 
remains is the electron's ability to undergo transitions 
between two discrete energy states. In order to fully un- 
fold the connection between the QED and the theory of 
two-level systems, we shall perform the second quantiza- 
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tion of the standard theory of qubits. The description of 
two-level systems in terms of creation and annihilation 
operators has been introduced before (of., for example, 
[3]) but no one has exploited the full potential of this 
formulation. The crucial new element in our formulation 
is the systematic use of Feynman diagrams. To expose a 
close analogy with the relativistic theory, including the 
form of the propagators, we shall choose the energy scale 
in such a way that the energy levels of the two-level sys- 
tem have opposite signs. In this way, we arrive at a pic- 
ture of a two-level system that coincides with the Dirac- 
sea view of quantum electrodynamics. The ground state 
of the two-level system corresponds to the occupation of 
the negative energy state, while the excited state cor- 
responds to the occupation of the positive energy state 
accompanied by a hole in the negative energy sea. The 
transition between these two states due to the interaction 
with a photon can be represented by the two elementary 
Feynman diagrams shown in Fig. [TJ 

There are significant advantages in using the Feynman 
diagrams and the Feynman propagators associated with 
these diagrams as compared to the standard perturbation 
theory used in nonrelativistic quantum mechanics. 

First, we never need the formula for the ground state 
expressed in terms of the noninteracting particles. This 
is due to the stability of the ground state under the 
adiabatic switching-on of the interactions. In the Feyn- 
man approach the difference between the physical ground 
state of interacting particles and the ground state of 
noninteracting particles amounts only to the phase fac- 
tor corresponding to all disconnected vacuum diagrams 

M. 

Second, a single Feynman amplitude combines several 
terms of the standard perturbation theory since in the 
Feynman approach all processes that differ only in the 
time ordering of the vertices are described by one Feyn- 
man amplitude (Fig. [2]). The number of diagrams of the 
standard perturbation theory that are combined into one 
Feynman diagram grows exponentially with the number 
of vertices. 

Third, there are many sophisticated tools available to 
evaluate and analyze Feynman propagators that greatly 
simplify the calculations and also give a deeper insight 
into the physical processes described by these propaga- 
tors. In particular, we shall use the quantum linear re- 
sponse theory to calculate the atomic polarizability and 
the spin susceptibility from the Feynman propagators. 
Our formalism is not restricted to two-state systems. It 
can easily be generalized to many-state systems (qudits) 
and we analyze as an example a four-state system — the 
atomic dipole — to show that the whole framework can 
easily be extended to cover this case. The main mes- 
sage of our investigation is that the Feynman description 
of quantum phenomena, known for its elegance, versatil- 
ity, and effectiveness in relativistic quantum field theory, 
also leads to significant simplifications in the theory of 
qubits. Of course, we are not trying to imply that qubits 
are relativistic objects. We shall only exploit formal sim- 



ilarities and use many available tools of a relativistic the- 
ory. Feynman propagators and Feynman diagrams in our 
approach should be treated as purely mathematical con- 
structs introduced as a means to streamline and organize 
perturbation theory. They greatly simplify the calcula- 
tions but they do not represent any physical objects. 

There is a huge number of papers and even a mono- 
graph [6j dealing with the theory of two-level systems 
and its applications. We believe that the point of view 
described in this paper will further our understanding of 
these systems. Our research has been prompted by a re- 
cent calculation of the atomic polarizability by Loudon 
and Barnett Q. Our results differ from their results in 
the fourth order of perturbation theory because they have 
not taken into account all the necessary corrections. The 
crossing symmetry of the polarizability, that played an 
important role in the derivation of the final result by 
Loudon and Barnett, is automatically satisfied in our for- 
mulation. In quantum field theory the crossing relations 
follow from the analytic properties of the propagators as 
functions of the energy parameter and from the direct 
connection between the polarizability and the retarded 
photon propagator. This connection enabled us to eas- 
ily calculate the polarizability of a two-level atom and 
the spin susceptibility in the fourth order of perturba- 
tion theory by evaluating the contributions from only a 
few Feynman diagrams. 

Our results clarify certain issues, like the opposite sign 
versus equal sign prescription or the darnping in the 
ground state, that are still being debated [1, H, [l^, |Tl|, 
[l2l . [isl [T^ . We show that both sign prescriptions are cor- 
rect but they apply to different physical situations. The 
equal sign prescription is appropriate for the scattering 
situation when we control the initial and the final pho- 
ton states. The opposite sign prescription is appropriate 
in the linear response theory when we control the initial 
state and also the form of the perturbation but we per- 
form a summation over all final states. Thus, only the 
opposite-sign convention is appropriate for the calcula- 
tion of the atomic polarizability. We also show that even 
though, as stated in |6|], "A two-level atom is conceptu- 
ally the same kind of object as a spin-one-half particle in 
a magnetic field" , the dynamical properties of these sys- 
tems are quite different. The differences become signifi- 
cantly different in higher orders of perturbation theory. 

Of course, one should keep in mind that our calcula- 
tions of atomic polarizabilities should not be taken too 
seriously because the two-level model gives only a very 
crude description of a real atom. However, for a single 
spin system, our results are close to reality. The only 
approximation being made in this case is that the posi- 
tion of the spin is frozen — the translational degrees of 
freedom are suppressed. 

It has been fully recognized that quantum field the- 
ory would, in principle, give unambiguous answers to all 
such questions but the prevailing opinion that "there are 
considerable difficulties associated with the treatment of 
optical damping in a non-phenomenological manner" Q 
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FIG. 2: Two time orderings in the standard perturbation 
theory that are combined into one Feynman amplitude repre- 
sented by one Feynman diagram. 



discouraged efforts to apply field-theoretic methods. In 
this paper we show how to overcome these "considerable 
difficulties". We formulate a theory that is simple be- 
cause it follows all the rules of a well established theory 
and it also has an unambiguous interpretation because it 
is systematically derived from first principles. 

In what follows we shall use most of the time a conve- 
nient system in units in which h ~ 1 , c = 1 , and /ip = 1 • 
Of course, in this system of units also eo = 1- More pre- 
cisely, we express every physical quantity in powers of 
the meter and fi, c, /zq (or eg) and then we drop fi = 1, 
c = 1, /io = 1 and eo in the formulas. For example, 
the Bohr magneton in these units is ^xb — 5.847 10^^'' m, 
Tesla is IT = 5.017 10^^ m'^, and the elcctronvolt is leV 
= 5.068 lO^m-i. 



II. THE MODEL HAMILTONIAN 



The physical system that we shall have in mind is pri- 
marily a spinning electron trapped in a spherically sym- 
metric potential subjected to a constant magnetic field 
and interacting with the quantized electromagnetic field 
and possibly an external time-varying electromagnetic 
field. We find it convenient to call this system the elec- 
tron to stress the analogy with quantum electrodynam- 
ics although it is a highly reduced model of an electron. 
We shall treat in detail the spin system coupled to the 
electromagnetic field through its magnetic dipole but we 
shall also extend our analysis to atoms coupled through 
their electric dipole moments. There are two cases here 
that must be distinguished: the literal two-level atom 
that requires a two-dimensional Hilbert space and an 
atom with a true electric dipole moment that requires 
a four-dimensional Hilbert space that can accommodate 
the three-dimensional dipole vector. 

The Hamiltonian H = Hq + Hj for the spin system in 



the second-quantized form is 

Ho^ J d^ri>\r)H^ip{r) 

+ \jd\: (E^r) + B^r)) : , (la) 

Hi = -fi Jd^rip^{r)eTil>{r)-B{r), (lb) 

where is the quantum-mechanical Hamiltonian of the 
electron in the absence of the magnetic coupling and the 
colons, as usual, denote the normal ordering. We shall 
assume that the magnetic moment of the electron is cou- 
pled to a constant external magnetic field and to the 
quantized magnetic field. Next, we assume that only the 
spin degree of freedom is active. Therefore, we can re- 
tain only one term in the expansion of the electron field 
operator 

xP{r)^x{r)ip, (2) 

where x('") is a fixed orbital electron wave function as- 
sumed to be spherically symmetric. The two-component 
fermionic operators are t/?^ = (V'liV'J) f^^id ip = {ipf,,tpg)- 
Their components create and annihilate the electron 
in the upper (excited) or lower (ground) energy state. 
Within this approximation, the Hamiltonian can be 
rewritten in the form 

Ho = fiBoiP^a^xP + i Jd\ : {E^{r) + B\r)) : , (3a) 

Hi = -^il)^(Txl}-Jd^rp{r)B{r). (3b) 

The parameter /i is the magnetic moment, Bq is the 
constant magnetic field (pointing in the z-direction) , 
and (T = (cj;, , (Tz) are the three Pauli matrices. In 
the interaction Hamiltonian the magnetic field opera- 
tor B is averaged with the electron distribution function 
pI*") = X*('')x('') over the region where the trapped elec- 
tron is localized. 

The Hamiltonian H = Hq + Hi conserves the num- 
ber of electrons. It acts independently in each subspace 
with a given number of electrons. Since there are just 
two creation operators in this model, the electronic Fock 
space is four dimensional. It comprises a one-dimensional 
zero-particle subspace, a one-dimensional two-particle 
subspace, and a two-dimensional one-particle subspace 
spanned by the state vectors V'I|0) and ipl\0). This two- 
dimensional subspace will be our qubit space. The stan- 
dard fermionic anticommutation relations 

{^„^,} = 0, {Vl,^]}=0 (4) 

imply that the operators i/j^cti/j annihilate the zero- 
particle and two-particle sectors, whereas in the qubit 
space they act as the Pauli matrices. Therefore, in the 
qubit subspace the Hamiltonian ([3]) is equivalent to the 
following one obtained from ([3]) by replacing all bilinear 
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combinations xjj'^aixjj of the operators ■0^ and xjj by the 
corresponding Pauli matrices: 

Ho = -fiBoa, + ljd^r: {E^r) + B^{r)) : , (5a) 

Hi = -f^o-J (frp{r)B{r). (5b) 

To stress the analogy between QED and quantum elec- 
trodynamics of two-level systems, from now on we shall 
denote the energy (iBq by the letter m. 

A. Spin system as a dimensional reduction of QED 

The formulation that employs the electronic creation 
and annihilation operators will enable us to define new 
objects — the propagators — that do not appear in the 
standard description of a spin system. The electron prop- 
agators, being auxiliary objects without direct physical 
interpretation, fully deserve the name "dead wood" , as 
Dirac 15] called them. However, a complete formulation 
of QED (including renormalization) without the prop- 
agators would be extremely complicated, if possible at 
all. We shall show that they are also very useful in the 
description of two-level systems. 

The Hamiltonian ([3]) acts independently in each sec- 
tor with a given number of electrons, but the electron 
creation and annihilation operators cause transitions be- 
tween these sectors. This leads here, like in full QED, to a 
greater flexibility of the mathematical formalism and will 
allow us to introduce objects that are not available in the 
standard theory of qubits based on the Hamiltonian ^ . 
Long time ago, the same idea has been successfully ap- 
plied to the study of the Ising chain [l^ and that served 
as an inspiration for the present research. The represen- 
tation of the spin operators as bilinear expressions of the 
creation and annihilation operators is the key ingredient 
of our approach. It enabled us to introduce the fermionic 
Feynman propagators and to employ the Wick theorem 
in its most convenient, field-theoretic form that leads di- 
rectly to standard Feynman diagrams. In contrast, the 
use of the spin operators as basic variables, does not lead 
to the Feynman rules in their simplest form known from 
QED. 

In order to better explain the relation between QED 
and our treatment of two-level systems, let us observe 
that the Hamiltonian ^ can be obtained by the dimen- 
sional reduction from three to zero spatial dimensions. 
To carry out this reduction, we drop entirely the coordi- 
nate dependence and we disregard the integration in the 
QED Hamiltonian Hd of the Dirac field 

Hd= j (fr {c'il)\r)csL-p-ilj{r) +rac^-ip\r)l3i}{r)) . (6) 

We keep only the mass term and we replace the Dirac 
field operator (V'i(r), ■02(t"), V'3('')j ■'/'4('')) by the space- 
independent operators {ipeT'4'g)- The operator i/)e anni- 
hilates the particle in the positive energy state and ipg 



annihilates the particle in the negative energy state. The 
rest energy toqc^ of the electron is to be identified with 
^Bq- Despite these drastic simplifications, we shall still 
retain the full analogy with quantum electrodynamics. 
This will enable us to use the highly developed formal- 
ism of QED and also to gain deeper insights that go with 
it. 



B. Magnetic dipole Hamiltonian 

Under the assumption that only the spin degree of free- 
dom is active and the orbital part of the electron wave 
function x(r) is fixed and spherically symmetric, only the 
magnetic dipole component of the radiation field is cou- 
pled to the electron. Therefore, it is most convenient to 
employ the multipole expansion, i.e. the decomposition 
of the electromagnetic field into the eigenstates of the an- 
gular momentum. Then, the integration of the magnetic 
field vector with the spherically symmetric distribution 
in the interaction Hamiltonian (j3bp eliminates all multi- 
poles except the magnetic dipole. We present the details 
of this calculation in the Appendix \X[ We shall rewrite 
the Hamiltonian (|A13p derived there as follows 

/•OO 

H = mij}^ a zip + I dk uj c\{k)ci{k) 
, Jo 

poo 

+ i>''(Tip- dkg{k)cf){k), (7) 
Jo 

where we introduced the dipole vector field 4){k) built 
from the Cartesian components of the annihilation and 
creation operators 



c^{k) + cl{k) 



2k 



(8) 



The formfactor g{k) is defined in Eq. (|A16|) and according 
to the formula (jA17[) it is proportional to the Fourier 
transform p(k) of the distribution function p{r) 



p{k). 



(9) 



The normalization condition imposed on p requires that 
p(0) = 1. Therefore, for small values of k the formfactor 
behaves as g{k) « pk^ j'K^fZ. To illustrate this property, 
let us consider the qubit realized as the spin degree of 
freedom of a nonrelativistic electron in the ground state 
of the Coulomb potential. In this case the distribution 
function p(r) and the corresponding formfactor .g(fc) are 



1 

pk^ 



1 



7rV3(l + fc2a2/4)2' 
where ag is the Bohr radius. 



(10a) 
(10b) 
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The applicability of the model interaction Hamiltonian 
d?]) extends beyond the simplest case considered here. 
Should the distribution function p(r) be of a more general 
character or the internal degrees be more complicated, 
the elimination of higher multipoles could still be justified 
as an approximation based on the small value of the ratio: 
atomic size/ wavelength. 



C. Two-level atom Hamiltonian 

In the case of a literal two-level atom considered by 
most authors, only one component of the electromagnetic 
field is coupled to the atom. Namely, the component that 
causes transitions between the ground state and one se- 
lected excited state. Therefore, it is sufficient to replace 
the three-component vector 4>{k) by a single component 
(j){k). In this way, we obtain the standard Hamiltonian 
for a two-level atom interacting with the quantized elec- 
tromagnetic field in the form [j, 01 

H — TTUTz+'S^ / dk Lo c\(k)ci(k) 

poo 

+ / dkg{k)(j,{k), (11) 
Jo 

which after the second quantization becomes 

/•OC 

H = mip^ azip + 'y^ / dk LO c\{k)ci{k) 

. "'0 



+ / dkg{k)(f>{k), 



where the formfactor g{k) 



m 



djf_ 



h{k). 



(12) 



(13) 



is obtained from the formula ^ by replacing the mag- 
netic dipole PL and its distribution function p by the elec- 
tric dipole d and its distribution function k. This natu- 
ral prescription will be confirmed in the next subsection 
when we derive the interaction Hamiltonian for a true 
atomic dipole vector. We place a hat on the symbols of 
all quantities that refer specifically to two-level atoms to 
distinguish them from the corresponding quantities for 
the spin system. 



three states corresponding to the degenerate upper level. 
This is precisely the situation in real atoms if the tran- 
sitions take place between the ground S state and the 
three excited P states. The inclusion of all three P states 
leads to full rotational invariance. Using this specific ex- 
ample we show how to extend our formalism to A^-level 
systems. The Hamiltonian H — Hq + Hj expressed in 
the formalism of second quantization can now be written 
in the form (cf. Appendix [X| 

/"OO 

= t/^tmt/) + V / dkwd\{k)di{k) 



+ i>Wip- dkg{k)4>{k), 
Jo 



(14) 



where we kept the same symbol ^(fc) to denote the elec- 
tromagnetic field because the change from the magnetic 
dipole field to the electric dipole field does not change 
any of the mathematical properties of the field cf){k). We 
introduced four annihilation and four creation operators 
corresponding to four atomic states. The operators for 
the ground state and the operators for the excited states 
in the Cartesian basis are combined into four-dimensional 
objects ip = {^^,^y,i;z,iJg} and ip'^ = {V;!, V'l, 
They obey the fermionic anticommutation relations ([4]). 
The matrices m and t are defined in Eqs. (|A22[) . The 
derivation in Appendix [Al of the formula for the formfac- 
tor function g(k) gives the precise meaning to the dipole 
moment d of the atomic transition and the dipole distri- 
bution function K(r) and its transform k{k). 



m 



de 



K{k). 



(15) 



Since for small values of k we have ji(fcr) « kr/i, the 
function K(fc) has the same normalization as p{k) — it 
approaches 1, when /c — > 0. In particular, for the P-S 
transitions in the hydrogen atom we obtain 



K{r) 



A'Ka^d^/& 
dp 



exp 



3r 
2ao 



7r\/3(l + 4Pa2/9)^' 
2i5/2eao 
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(16a) 
(16b) 
(16c) 



D. Electric dipole Hamiltonian 

The truncation of the atomic Hilbcrt space to only 
two dimensions does not allow for the construction of an 
atomic dipole vector that could be coupled to the electric 
dipole field. Such a construction can be carried out if we 
enlarge the Hilbert space of the relevant atomic states 
to four dimensions. We shall still have only two energy 
levels but in addition to the ground state we introduce 



E. Conservation of angular momentum 

The interaction Hamiltonian for the spin system is in- 
variant under all rotations since it is a scalar products 
of two vectors. However, the full Hamiltonian is invari- 
ant only under rotations around the z axis since the free 
fermion Hamiltonian ([7]) contains the z component of the 
vector (T. The physical origin of the symmetry breaking 
is the external magnetic field fixed along the z-axis. It 
splits the energy levels of the magnetic dipole and breaks 
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the full rotational invariance. In contrast, the Hamil- 
tonian for the electric dipole is invariant under the full 
rotation group. This invariance is possible because the 
Coulomb potential of the hydrogenic atom is rotationally 
symmetric and wc have included all three components of 
the excited P state. These components form a vector 
representation of the rotation group. 

The invariance of the Hamiltonian implies the comniu- 
tativity of the angular momentum operator Mz with H 
leading to the conservation of the in both cases. The 
angular momentum operators for the spin system and for 
the electric dipole are 

1 f°° 
M,^ -iP^'a.ip-i J dke^:jkc]{k)ck{k), (17) 

M, = ^\k)ilj{k) - I j dke,jkd]{k)dk{k), (18) 

where the spin-one matrices Si with elements {si)jk — 
—ieijk act in the subspace of excited states. Conservation 
of angular momentum during interaction becomes obvi- 
ous when the angular momentum operator and interac- 
tion Hamiltonian are written in the angular momentum 
basis. We shall use the spin system to illustrate these 
properties. Let us construct the components of the mag- 
netic dipole field 4>±{k) and 0o(fc) from the annihilation 
and creation operators of photons with the definite angu- 
lar momentum Mz = ±1,0 introduced in the Appendix 



c_(fc) -cUfc) 



2k 



cUk)~c+{k) . 
'(k) = 7= = 0+(fc), 



2k 

h 



Mk) 



co(fc) +cT(fc) 



2k 



(19a) 
(19b) 
(19c) 



The operators Mz and Hj take now the form 



Mz = iv-V.V + / dk 4(fc)c+(fc) - cUk)c^{k) 



(20) 



Hi ^i^^'a+i: dkg{k)(l3^{k)+i;''a-^ dk g{k)4>+{k) 
Jo Jq 



V'VzV' / dkg{k)(f>o{k), 
Jo 



where 



V2 



V2 



(21) 



(22) 



The field <j)+{k) coupled to (t_ annihilates the photon 
with Mz = — 1 or creates the photon with Mz = 1. Thus, 
it increases the angular momentum by one unit. The field 
(j)- (fc) coupled to (T+ decreases the angular momentum by 
one unit. Each term in the Hamiltonian (j2ip conserves 



angular momentum. For example, when (t+ transfers the 
electron from the ground state to the excited state in- 
creasing its angular momentum by one (the first term), 
the angular momentum of the electromagnetic field de- 
creases by one unit. Similar analysis can be carried out 
for the electric dipole. Of course, for the literal two-level 
atom there is no invariance under rotation because only 
one angular momentum state of the photon interacts with 
the atom. Hence, only one component of the electronic 
P state (and not all three) can be excited. 



F. Time-reversal invariance 

Both theories, describing the spin and the two-level 
atom, are invariant under the time reversal. This invari- 
ance can be proven directly but it also follows from the 
fact that our models are obtained by the dimensional 
reduction from QED which has this property. Time- 
reversal invariance is an important requirement to obtain 
a correct description of the optical damping, as stressed 
in Ref. [10]. In what follows we shall make use of this 
invariance. Under the time reversal the signs of the fre- 
quency and angular momentum are reversed. Therefore, 
there is no need to calculate the photon propagator for 
the negative values of Mz for the spin system because 
they can be obtained from those for the positive values 
by reversing the sign of the frequency. When the results 
are the same for positive and negative values of Mz, as is 
the case for the atomic system, time-reversal invariance 
means that the photon propagator is an even function 
of the frequency. The conservation of angular momen- 
tum and time-reversal invariance simplify the calcula- 
tions since they reduce the number of Feynman integrals 
that are to be evaluated. 



III. PROPAGATORS AND THE S MATRIX 

All transition amplitudes can be expressed in terms of 
Feynman propagators — the expectation values in the 
ground state of the time-ordered products of field op- 
erators. Since we shall be working within perturbation 
theory, the most useful representation of the propagators 
is the one that is based on the perturbative expansion of 
the S matrix. The relevant formula for the S matrix is 
the following standard expansion into the time-ordered 
products of the interaction Hamiltonians [l3| : 



= Texp (^-i J dtHi{t) 



E 

n=0 



dh--- dtnT[Hi{ti)---Hi{tr,)]. 



(23) 



The interaction Hamiltonian in this formula is taken in 
the Dirac picture. We shall introduce all the necessary 
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theoretical tools starting with the spin system but later 
extending them to atoms by making obvious modifica- 
tions. We will find it expedient, even though it is not nec- 
essary since there are no infinities, to perform the mass 
renormalization. This amounts, exactly like in QED, to 
adding the mass-correction term Smip^azip to the free 
Hamiltonian and subtracting the same term from the in- 
teraction Hamiltonian. In our case, the freedom of choos- 
ing Sm can be viewed as a mechanism to improve the con- 
vergence of perturbation theory. After the mass renor- 
malization, the free Hamiltonian and interaction Hamil- 
tonian in the Dirac picture become 



dk:{TT^{k) + k^<l>^{k)) 



(24) 



Hi{t) ^ e'^"*iJ/e-*-f^°* 

= i)^t)(Tip{t)- dkg(k)cj}{k,t) - 5mip'' {t)azXp{t), 
Jo 

(25) 

where 7r(fc) is the canonically conjugate momentum 

<k)^~z^{cik)~cHk)). (26) 

The time dependence of the operators ip{t),'ip'^ (t), and 
(p{k, t) is determined by the renormalized fermionic 
Hamiltonian (|24p and it has the following form: 



'g 



(27) 
(28) 



where m = ttiq + Sm. The time dependence of the field 
4>{k,t) is 



</.(fc,t) 



c(fc)e-*"* -hct(/c)e''^* 



2k 



(29) 



Note that due to our normalization, the electromagnetic 
field operators 4>{k,t) and Tv{k',t) = (p{k',t) satisfy the 
equal-time canonical commutation relations 



[4>i{k, t), TTjik' , t)] — i6ijd{k — k'). 



(30) 



In order to describe the interacting system, we need 
the propagators defined in terms of the field operators 
^(t), ^''{t), and $(fc, t) evaluated in the Heisenberg pic- 
ture. We shall use lower case and upper case letters to 
keep the distinction between the Dirac (interaction) pic- 
ture and the Heisenberg picture operators. The Heisen- 
berg picture operators obey the following equations of 
motion: 



{idt - mocr^)*(t) 



dkg{k)(T-^{k,t)'^{t), (31a) 



{d^ + k^)^{k, t) = -g{k)-^Ht)cT-^{t). (31b) 

The canonical equal-time commutation relations of the 
Heisenberg operators are the same as their free counter- 
parts 



[^^{t)y^{t)]=5^p, (32a) 
$,(fc,t),$j(fc',t)j =i5^j5{k-k'). (32b) 



All remaining commutators or anticommutators vanish. 

The perturbation expansion of the propagators can be 
obtained from the following formula 0, [1^1 by expanding 
the time-ordered exponential function into a power series 
according to Eq. 



{G\T[^{t,) ■ ■ ■ ^{U)^\t\) ■ ■ ■ ^\t',)Hk,,t'i) ■ ■ ■ m,ti)\ \G) 

_ {g\T [V(ti) • • • ^UWiA) ■ ■ ■ i:Ht'Mki,t'{) ■ 



■<l)iki,t'{)exp{-iJdtHiit))] \g) 



{g\Texp{-iJdtHiit)) \g) 
I 



(33) 



We have omitted here all indices leaving only the depen- 
dence on time and on the wave vector. The operators 
on the left hand side of this equation are in the Heisen- 
berg picture while those on the right hand side are all in 
the Dirac picture. In this formula \G) denotes the true 



ground state of the interacting system and \g) denotes the 
ground state of the free Hamiltonian Hq. In the state \g) 
there are no photons and the negative energy state of the 
electron is occupied. The advantage of using this funda- 
mental result, already mentioned in the Introduction, is 
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that the detailed knowledge of the ground state \G) is 
not needed. The difference between the state vectors \G) 
and \g) is just a phase factor and the denominator in 
the formula (j33p representing the contributions from all 
disconnected vacuum diagrams takes care of that. 



IV. FEYNMAN DIAGRAMS AND FEYNMAN 
RULES 



In order to derive the Feynman rules that connect the 
Feynman diagrams with the corresponding transition am- 
plitudes we start, as in QED, from the free field opera- 
tors. The time evolution of these operators is given by 
Eqs. (P7| and (PI) . 

The basic ingredients of the Feynman formulation of 
QED are the free one-electron propagator Sp and one- 
photon propagator Dp. In our model they are defined as 
follows: 

SFo.0it - t') = -i{g\T (^^o^mlit')) \g), (34) 
Df,, (fc, k\ t-t') = -i{g\T (Mk, t)c^o{k\ t')) \g), (35) 



where Sf{po) has the form 



Sf{po 



Po-me + le Po 



POCTz 



m 
1 



Po - {m- ie)az ' 



(38a) 
(38b) 
(38c) 



The formula (|38ap holds also for the atomic dipole when 
the excited states form a subspace. In what follows we 
shall use the same symbols Pe and Fg to denote the pro- 
jectors in all three cases. It will be clear from the con- 
text, whether Pe projects on the one-dimensional sub- 
space (spin and two-level atom) or on the three dimen- 
sional subspace (atomic dipole). As compared with the 
Fourier transform of the electron propagator in the rela- 
tivistic theory 1/ ■ p— m + ie), the two-level propagator 
([55)1 lacks the spatial part of the momentum vector and 
has the Pauli (Tz matrix instead of 70. The presence of 
az in the numerator in Eq. (|38ap reflects the fact that 
we work with ■0t instead oi — "070. We shall use 
the same symbols to denote the propagators and their 
Fourier transforms. The arguments will always indicate 
which is the case. 



where \g) is the ground state of the system without inter- 
action. We have introduced the photon propagator only 
for those photons that are coupled to the electron. 



A. Free electron propagators 

The free electron propagator is easily evaluated with 
the use of Eqs. ^7} taking into account that the only 
nonvanishing matrix elements of the bilinear product of 
the creation and annihilation operators are {g\'4'e'4'l\9) = 
1 and {gl'ifj'j^ipglg) = 1. Therefore, we obtain 

iSFap{t - t') 

=e{t - t'){g\Mmlit')\9)-0it'- t){g\i;l{t')Mt)\9) 



--e{t - t')Pea/3e-""=(*-*') - e{t' - t)Pgc,/3e-™«(*-*'\ 



(36) 



where Pe = (1 + az)/2 and Pg = (1 - (Tz)/2 are the 
projection matrices on the upper and lower energy states, 
respectively. For the spin system and the two-level atom 
we have nie — m and rUg = —m. However, for the dipole 
atom these two parameters will be independent. The 
final result can be expressed in matrix notation (omitting 
the indices a and (3) as the following Fourier integral: 



SF{t-t') 



dPO a I \ ~ 

— SF[Po)e 
, 2tt 



-ipo (*-*') 



(37) 



B. Free photon propagators 

The free photon propagator is 

DF^J{k, k\t- t') - -te{t - t'){g\Uk,mj{k\t')\g) 
-i0it' ~t){g\Uk',t')^j{k,t)\g) 

■^^Ak^^^-^uJ\t-t'\ 



-'>■— — tt; -e 

2k 



We shall also need its Fourier representation 



(39) 



DF^J{k,k',t-t') = / '^^DF^,{k,k\ko)e-'^"'^'-''\ 



where 



DFr,{k,k',ko)= 
d^j5{k - k') 



2k 



(40) 

(41a) 
(41b) 



All Feynman amplitudes can be constructed from the 
electron propagator ((37| . the photon propagator (|4ip . 
the vertex, and the mass insertion following the same 
general rules as in QED. The starting point is the defi- 
nition p3p of a general propagator. In the n-th order of 
perturbation theory the contribution to the propagator 
is expressed as an expectation value of the time-ordered 
product of operators ^p, and <j) integrated over n time 
variables. In our model, as in the standard QED, all these 
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expectation values can be evaluated with the help of the 
Wick theorem (cf., for example, The only dif- 

ference in applying this theorem is, in contrast to QED, 
that we have not interchanged the creation and annihi- 
lation operators for the negative energy state. Calling 
the electron in the ground state an antiparticle would 
stretch the analogy with QED too far. Therefore, in our 
case the normal ordering means that all operators ipl and 
4>g stand to the left of all operators -0^ and ipe- 



C. Feynman rules 

The scattering amplitudes in QED are commonly eval- 
uated in momentum representation. In our case, the 
transformation to momentum representation means the 
transformation from the time domain to the frequency 
domain. The Feynman rules in the frequency domain 
are obtained by substituting everywhere the free electron 
propagators and photon propagators in the form of the 
Fourier integrals (|57|) and (|¥T|) . Next, in the n-th order 
of perturbation theory we perform n time integrations. 
Finally, we take the inverse Fourier transforms with re- 
spect to all remaining time arguments of the propagator 
([55)1 . These operations lead to the following Feynman 
rules: 

• Each electron line corresponds to the Fourier trans- 
form of the electron propagator and is represented 
by iSpipo)- 

• Each photon line corresponds to the Fourier trans- 
form of the photon propagator and is represented 
by iDptjik, k' , fco). 

• Each vertex is depicted by two electron lines and 
the photon line meeting at one point. It is repre- 
sented by —iVi{k) — —ig{k)ai. The energy conser- 
vation at each vertex results in the appearance of 
2Tr5{po - go - ^o)- 

• Each mass insertion is depicted by a cross where 
two electron lines meet. It is represented by iSmaz- 
The energy conservation at each mass insertion re- 
sults in the appearance of 2Tr5{po — qo)- 

• All 2x2 matrices corresponding to electron propa- 
gators are multiplied in the order indicated by the 
arrows on the diagram. 

• Each closed electronic loop brings in a minus sign 
and a trace over the matrix indices. 

• There is a summation over all repeated vector in- 
dices and an integration over all repeated values of 
the length of the wave vector. 

• There is one integration over the energy variable 
for each closed loop, accompanied by the division 
by 27r. 



po 




iSp (po)a/3 
iDpijik, fc', fco) 



-iVi{k) 



af3 



po 



-X- 



FIG. 3: Feynman rules. For clarity, we have written explicitly 
all indices. 



These rules are summarized in Fig. [3l Calculations of 
the lowest order radiative corrections to the electron and 
photon propagators based on these rules are presented in 

Sections EH and IVnl 

In the case of the two-level atom the only changes in 
the Feynman rules as compared to the case of the spin 
system is that the elementary vertex is represented just 
by —ig{k)<Tx and the photon propagator has no indices. 
In the case of the atom with an electric dipole the free 
photon propagator retains its form ([1T|) . The free elec- 
tron propagator must be taken in the general form (j38a[) 



(42) 



and at each vertex the matrices cr must be replaced by 
the matrices r. 



V. RADIATIVE CORRECTIONS 

Owing to the absence of the space components of mo- 
mentum vectors, the calculation of radiative corrections 
is much simpler here than in the full-fledged QED. There 
is no need to combine denominators a la Feynman and 
Schwinger. All integrations with respect to the loop vari- 
ables Po,ko etc. can be evaluated analytically by the 
residue method in any order of perturbation theory. At 
the end we will be left only with the integrals over the 
wave vectors of photons weighted with g^(fc). Of course, 
those integrals cannot be evaluated if the function g{k) 
is not specified. 

In order to explain how the calculations are done, let 
us consider an integral represented by an arbitrary Feyn- 
man diagram. The integrand is a product of electron and 
photon propagators. To perform all the integrations with 
respect to the loop variables, one may choose the electron 
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(a) 




(c) 



— ^ 



(d) 



(e) 




FIG. 4: Feynman diagrams representing the lowest-order ra- 
diative corrections to the electron propagator, photon propa- 
gator, and the vertex part. 



(a) 



(b) 









— ^— 





FIG. 5: Graphical representation of the relationship between 
the propagators and the corresponding self-energy parts. The 
double-lines represent full propagators and the gray box and 
circle represent the self-energy parts. 



propagator in the form (|38ap and use the photon propa- 
gator in the form (I41bp . The numerator of the integrand 
corresponding to each Feynman diagram is a polynomial 
in the integration variables. The denominator is a prod- 
uct of first-order polynomials in the integration variables, 
each factor leading to a simple pole. All integrations can 
easily be done by the standard residue method. Note 
that after each successive integration the integrand re- 
tains its rational form. Therefore, it will continue to be 
amenable to the same treatment as during the first in- 
tegration. Alternatively, we may choose the interaction 
Hamiltonian in the angular momentum basis (|2ip . The 
following algebraic properties of the matrices a± are then 
very useful: 

crl=0 = , cr+CT_ = 2Pe, cr-cr+ = 2Fg, (43a) 
(t„Mct„ = (T+Ma- + a^M(T+ + (J^^Ma^,, (43b) 

n 

^a„PeC7„ = Pe + 2P3, ^<7„P3a„ = Pg + 2Pe, (43c) 
n n 

where M is an arbitrary matrix. With their help, and 
using the anticommutativity of a± with cr^ , we can reduce 
every Feynman integral to to a very simple form. 

In the case of a two-level atom the calculations are sim- 
pler than in the case of the spin system. Due to the ap- 
pearance of only the matrix in all vertices, the matrix 
algebra is almost trivial. In each integrand we can bring 



up front all ax matrices using the relations (TxCz = —o'zO'x 
and (7^ = 1. Therefore, each time we interchange the or- 
der of ax and az in the electron propagator the sign of 
az must be reversed. Since there will be an even number 
of vertices in all the diagrams under consideration, the 
matrices ax will disappear completely and we will be left 
with a diagonal matrix that contains only the matrices 
az- The trace of such an expression is the sum of the 
terms corresponding to the eigenvalues ±1 of az- 

In the case of the electric dipole, the following algebraic 
properties of the r matrices: 



+ 3Pg, (44) 



used in conjunction with the general form (j38ap of the 
free electron propagator greatly reduce the number of 
integrals that are to be evaluated. 

We shall show how these rules work in practice by cal- 
culating radiative corrections to the electron and photon 
propagators. The procedure employed very often in QED 
relates the full electron and photon propagators to the 
self-energy parts. This procedure enables one to go be- 
yond the simplest version of perturbation theory and sum 
up an infinite (geometric) series. The self-energy is the 
sum of contributions from strongly connected diagrams, 
i.e. the diagrams that cannot be disconnected by cutting 
only one line. The relations between the full propaga- 
tors and the self-energy parts are shown schematically in 
Fig. El 



VI. ELECTRON PROPAGATOR 

In the case of the electron propagator Gf{po) the re- 
lation between the propagator and the self-energy part 
E(po), iUustrated in Fig. [5^, reads 



Gf(po) = Sf{po) + 5f(po)S(po)Gf(po)- 



(45) 



All three objects that appear in this equation are 2x2 
matrices. The iterative solution of Eq. that shows 
explicitly the relation between the propagator and the 
self-energy part, is 



Gf{po) = Spipo) + S'f(po)S(po)S'f(po) 

+ 'S'f(po)S(po)S'f(po)S(po)'S'f(po) 



(46) 



This formal geometric series can be summed up to the 
following compact form: 



Gf(po) 



1 



(47) 



where the inverse is to be understood as the inverse of a 
matrix. The series (1461) without resummation is mean- 
ingless because it is divergent when po ~ ™- 

The radiative corrections to the electron propagator in 
the second order of perturbation theory are represented 
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by the three Feynman diagrams (a)-(c) shown in Fig.|4l dipole atom &'^\po), constructed according to the rules 
The self-energy parts in this order for the spin system stated in the previous section, have the form 
S'^^)(po), for the two-level atom S'^^)(po), and for the 



^T.J ^^Y.j dk%{k)SF{po + ko)V,(k')Dp,,{k,k\ko) 

/OO 1 poo poo 

/ rffcV / dk'i:r{V,{k)SF{po)}Vj{k')DF,j{k,k',Q)-5ma,, (48a) 
-OO ■ JQ ■ Jo 



3 

roc 77, /.OO poo 



E(2)(po) = E(2")(po) + E(2-)(po)=i / ^/ dk I dk'V(k)SF{po + ko)V{k')DF{k,k',ko)~Sma,, (48b) 



27r 

OO 



^0 



/OO 11 pOO pOO 

T^E / '^^E / dk%{k)SF{Po + ko)V,{k')DF.,{k,k',ko)~5m. 
-OO j Jo J Jo 

(48c) 

The tadpole diagram (Fig. [4Jd) does not contribute in the case of the two-level atom and the dipole atom because 
Tr{axSF{pQ)} = and Tr{TiSF{po)} = 0. The analytic expressions for the self-energy parts obtained by the 
application of the Feynman rules are 

S(^Hpo)=^E rf^^» ^7 ] r—a , I 

^ Jo J-ca 27r pa + ko- (m- leja^ k^ - k^ + le 

-z^ rdkg^k) r^Tr(a„ / . , U„ J . -dma., (49a) 

„ Jo J-oo 27r [ Po- (m ~ iqcJz } -k^ + «e 

/'C'O pOO 11, -1 -1 

-:7^T„ ^ — + ^, ' r„ -Jm, 49c 

^TT \po + fco — ?Tle + 26 PQ + ko — TTlg — 16 J k^ — k^ + 16 



t^^\po)^iY.j^dkg\k) j 



where 5rh is the mass renormalization matrix with the eigenvalues 5nie and 5mg . With the use of the relations 
and (|44p we can replace all matrices by the projectors 

E '^"^ 1 r-—an = _^ . + ——r (50a) 

^ Po + ko ~ [m — t£)az po + ko + m — le po + ko - m + le 

{'""po _ ( ie)a J = { 2m{pl - 4m2 + ze)-i ^ = z) ' (^^^) 



1 



' Po + ka — [m — ie)az po + ko + m — ie p^ + k^ — m + ie 



(50c) 



/ T« ; H ; Tn = ; H : (50d) 

^ \po + ko — nie + le po + ko — nig — te J po + ko — irig — le po + ko — nie + le 



and then we can easily perform the integrations over kg (to> will be equal either to or nig) 
dko I 1 1 . f°° dpo 2m 



27r Po + ^0 — iTT'X T — k"^ + ie 2k{po ± k — ni\ ^ ie) ' 7_oo ^tt Pq — Am? -\- it 
Finally, we obtain 



(51a) 



r 



where 



Note that the contribution proportional to m^, corre- 

(53) 
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spending to the tadpole diagram, has the same form as 
the contribution from the mass correction. 
For the two-level atom, we obtain 



°dk g^k) 
'o 2k pq + k + m — ie 
'dk g^{k) 



2k po — k — m + ie 



■ e g 



(54) 



)6m. (55) 



The electron self-energy part for the dipole atom is 
slightly more complicated 



t^^Hpo)=K 



^dk f{k) 
2k Pq + k — nig — ie 

dk f(k) 

a I TTi T r~ PeOm-e - ¥g6mg. 

Iq 2k Pq — k — rrie + le 



(56) 



(57) 

The mass corrections Sm, 5rh, and Srh will be chosen so 
that the propagator G{po) with radiative corrections has 
a pole at the renormalized mass. These pole conditions 
imply that l^^^^maz) ~ and &^^maz) — and they 
give 



Sm 



dk 
2k? 



2 s 3fc -I- 2m 



rdk^^.^. 1 

/ ^5 fc) , , o ■ 
Jq 2k k + 2m 



(58a) 
(58b) 



For the dipole atom the mass corrections are different for 
the ground state and for the excited state — the energy 
of the excited state is raised and the energy of the ground 
state, as is always the case, is pushed down 



SlTle 



dk g^{k) 



2k k + Am 



dk g\k) 

om„ ^ —6 — — , 

^ ' 2k k + Am 



(59a) 
(59b) 



where Am 



me 



All these mass corrections give 



frequency-independent shifts in the level separation. The 
electron propagators do not have a direct physical inter- 
pretation but they serve as important ingredients in the 
calculation of the photon propagators. In particular, we 
will need the mass corrections to complete the calculation 
of the spin susceptibility and the atomic polarizability in 
the fourth order of perturbation theory. 



VII. PHOTON PROPAGATOR 

The photon propagator plays a distinguished role in 
our formulation, much more so than the electron propa- 
gator, since it enables one to calculate several important 
physical characteristics of two-level systems. The propa- 
gation of photons is, of course, modified by the presence 
of a two-level system. The scattering of photons off a 



two-level system is the counterpart of an important phe- 
nomenon in QED — the vacuum polarization. 

The relation of the full photon propagator to the self- 
energy part Ilij{k,k',ko) is illustrated in Fig. [SJs. It is 
slightly more complicated than in the case of the elec- 
tron propagator because, in addition to a multiplication 
of matrices in the space of the vector components, we 
must perform an integration over the wave vector k. The 
counterpart of Eq. (|45)) is 

Gpijik, k\ ko) 

/•OO />oo 

= L'fij(A:, fc', fco) + ^ / dki^^ dk2 
I Jo „ Jo 

X Dpii {k,ki, ko)Ilin (fci , fc2, ko)gFnj {k2 ,k',ko). (60) 

Taking into account the fact that DFij{k,k',ko) is pro- 
portional to the Kronecker 5ij and the Dirac 6(k — k'), 
we can rewrite this equation in the form 



GFij{k,k',ko) 



S.jSjk - k') 

/cq — fc2 + ie fcg — fc2 -I- ie 



xV / dk" g{k")Va(ko)gFij{k",k',ko), (61) 

where we took advantage of the factorization of 
n.y(fc, /c',A:o) 

U,,{k, k\ ko) = g{k)V,j{ko)g{k'). (62) 

The iteration of Eq. ((6T|l leads to the following expansion: 

gFi][k,k' ,ko) 

_ 5,,5{k-k') , g{k) g{k') 



kl-k'^ + ie kl-k'^ + 



le 



kl - fc'2 + ie 



k^ — k'^ + ie 



'^0 



fc2 - k"^+ie 



X Vijiko) 



9{k') 



kl - fc'2 + ie 



(63) 



This geometric series can be summed up and the final 
formula is 

GFijik^k'.ko) ^ DFij{k,k',ko) 

1 9{k) rj. n ^ 9{k') 

+ l^r^WTJf^^^^'^l^l^k^e 

The transition matrix T{ko) has the following represen- 
tation in terms of the self-energy part: 



T{ko) 
where 



Viko) 



1 



l + VikoMko) P(fco)-i + /i(fco)' 
dkg^{k) 



hiko) = 



(65) 



(66) 
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Both T(fco) and 'P(fco) in Eq. ([65)) are to be treated as 
3x3 matrices and the matrix to the power of —1 is meant 
as the inverse matrix. 

The function h{ko) wiU play an important role in our 
calculations because in the lowest order of perturbation 
theory its real part determines the shift in the position 
of the resonance and the imaginary part determines the 
width of the resonance 



Re h{ka) 
Im h{ko) 



fcg 



^g^(lfcol) 

2|fcnl ■ 



(67) 
(68) 



It follows from the assumptions that determine the valid- 
ity of our model that the real part of h(ko) is practically 
constant and can be replaced by its value at and the 
imaginary part varies as fcg. For example, when p(r) and 
g{k) are given by Eqs. pU]) . we obtain 



Hko) = 



(1 + „ 9^4 _ ^6 _^ iQ^^S 



127rai^(l+^2)4 



(69) 



where ^ = kQao/2. The value of the dimensionless pa- 
rameter ^ is very small in the range of wave vectors that 
cause the transitions between the two energy levels of 
our qubit. Thus, we can take only the leading terms and 
neglect higher powers of ^ as compared to 1, to obtain 



Hko) 



+ i 



(70) 



The formulas ([Mjl and ([65]) are also valid for the two- 
level atom and the dipole atom. In both cases h{ko) is 
defined by Eq. ([M]) where g{k) should be replaced either 
by g{k) or by g{k). Of course, in the first case there are 
no vector indices — T{ko) and P(fco) are not matrices 
but ordinary functions. In the second case, as is seen in 
Eq. ([75)1 below, owing to the full rotational invariance, 
the matrix V{ko) is proportional to 5ij. 



A. Second order of perturbation theory 

In the second order, the radiative correction to the 
photon propagator is represented by the diagram (a) in 
Fig. El The photon self-energy part, constructed accord- 
ing to the rules given in Fig. [3] has the form 



Tr {(TqS'f(po + ko)abSF{pQ)} ■ 



(71) 



The indices (a, 6) may take the values x, y and z in the 
Cartesian basis or the values — and in the angular 
momentum basis. The matrices cr are to be replaced by 
ax for the two-level atom and by the matrices t in the 
case of the atomic dipole. 



For the spin system it is convenient to choose the in- 
teraction Hamiltonian in the angular momentum basis 
([2T1) because in this basis the photon self-energy part is 
diagonal. The components of the self-energy in the an- 
gular momentum basis are 7'±(fco) and Vo^ko). They 
correspond to the following choices of the matrices a in 
Eq. (ini): 



V+iko) : 
Vo{ko) 



On 



G 



+ 



CTa = cr+ , 0]y= 0-_ 
CTa — Oz, Uh = (Tz- 



(72a) 
(72b) 
(72c) 



Making use of the properties of the a matrices, we 
end up with the following integrals: 



v'i\ko) 



dpo 



1 



1 



27ri po -I- fco T w ± ie po i TO T *e 
2 



2m =p fco ' 



V^'\ko) = 0. 



(73a) 
(73b) 



The component (^o) vanishes because in the corre- 
sponding integrals both residues lie in the same half- 

'P^^\~ko) is a direct 



(2) 

plane. The relation V_ {ko) 



confirmation of the time-reversal invariance. The an- 
gular momentum components of the transition matrix 
T(fco) obtained by substituting these self-energy parts 
into Eq. 



are 



(fco) 
To^'^(fco) = 0. 



2rn T fco — 2h{ko) 



(74a) 
(74b) 



For the two-level atom we must take a = x and 6 = a; in 
Eq. ([TT]) . After evaluating the trace, the integral reduces 
to the sum of two simple integrals 



J-oo 27r 



X Tr 



1 



Po + ko- {m 
dpo 1 



Pq 
1 



(to - ie)az 



_^ 2ni po + ko + m - ie po 
dpo 



1 



TO 
1 



, 27ri pa + ko — m + ie Po 
The result of the integrations is 



V^'\ko) 



4to 



4to2 



fc2 



(75) 



(76) 



and it leads to the following formula for T(fco) in the 
lowest order of perturbation theory: 



4to 



4to2 — ko — Amh{ko) 



(77) 
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(a) 



(b) 





(d) 



(e) 




(g) 



(h) 




FIG. 6: Photon self-energy diagrams in the second and fourth 
order of perturbation theory. 



For the dipole atom the contribution represented by 
the diagram (a) in Fig.[S]leads to the following expression 
for the self-energy part: 



dpo 



1 



1 



' dpo 1 1 

o 27ri po + ko ~ nig — iepo — vie + *e 
2Arn 



This leads to the transition matrix of the form 

2Am 



Ato2 -kl~ 2/\mh{ka) 



(78) 



(79) 



fco ~ 2m fourth-order terms dominate over second-order 
terms. The remedy comes from the realization that these 
double poles simply indicate an additional shift in the po- 
sition of the resonance. Indeed, expanding \/[2m—k^—8) 
into powers of 5 we encounter higher-order poles. We en- 
countered the same problem in the expansion of the elec- 
tron and photon propagators into the perturbation series 
but the difference is that then we were able to sum up the 
whole geometric series. Here, we can do it only order by 
order. In the present case, we can eliminate the double 
pole in the formulas (jB8[) by the following substitution: 



1 



2m — ko (2m — /cq)^ 2m — fco — 5 ' 



(80) 



that reproduces correctly the lowest order correction in 
d. Higher powers of 6 contribute to higher orders of per- 
turbation theory. Applying this procedure to the expres- 
sions for the self-energy parts (jB8p we obtain the fol- 
lowing formulas that do not suffer from the double-pole 
problem: 



v^i^'\k„) 



(2+4) 




(fco)=n^(fco) 
'dk g^{k) 



2(1 -b) 
2m =p /co ~ (5 ' 

(4)^ 



1 



/o k k + 2m [k + 2m)2 — k^ — ie 
Therefore, the transition matrix in this order is 



T. 



(2+4) 



(ko) = - 



2(1-6) 



2m=Ffco-'5-2(l-5)/i(fco)' 



(2+4) 



(fco) 



_ -0(2+4) 



(ko) 



(81a) 



(81b) 



(82a) 
(82b) 



The last equation follows from the fact that, as seen from 
Eq. (|73bp . Vo does not contain terms of the second order. 

The results for the two-level atom in the fourth order 
are even simpler since there is only one component of the 
self-energy part and there are no double poles. Substi- 
tuting the self-energy part (|B13p into the formula (|65p . 
we obtain 



B. Fourth order of perturbation theory 

The calculation of the photon self-energy part up to 
the fourth-order of perturbation theory requires the eval- 
uation of all the contributions to the photon self-energy 
represented by the Feynman diagrams (b)~(h) shown in 
Fig. [51 These calculations are presented in Appendix 
IB] Upon substituting the results of this calculation into 
Eq. (j65p . we obtain the formula for the transition matrix. 
However, there is an additional problem now that has not 
been present in the calculation of the self-energy part 
in the lowest order. In the formulas for the self-energy 
parts (jB8ap of the spin system we encounter double poles 
l/(2m — fco)^ and l/(2m -I- fco)^. Such terms indicate a 
breakdown of the simple perturbation theory since for 



f^'+'\ko) 



4m(l - b) 



Am^ 



ftp 



4m(l - b)h{ko) 



(83) 



The transition matrix for the dipole atom in the fourth- 
order, obtained from the self-energy part (jBiep . has the 
same general form as that for the two-level atom 



Tt''iko) = -S., 



2Am(l - b) 



Am2 - kl - 2Am(l - b)h{kn) 



(84) 



VIII. PHOTON SCATTERING AMPLITUDE 

The photon scattering amplitude (w) can be ob- 
tained [3i from the photon propagator ([M]) by strip- 
ping off the free propagators at both ends and putting 
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the whole expression on the energy shell kg = cu, k = 
uj, k' = oj. Therefore the scattering amplitude is related 
to the transition matrix by the formula 



(85) 



The argument w of the scattering amplitude can only take 
positive values because the photon energy is positive. 

In the second order of perturbation theory, the self- 
energy part for the spin system is given by Eq. ([73|) . 
Therefore, according to Eq. ([85]). the photon scattering 
amplitude for a spin system is 



2m — A(w) ^ Lo — iT^Lu) ' 



(86a) 
(86b) 



where the energy-dependent shift and width in this order 
according to Eqs. (IST)) are 



A(a;) = P 



'dk 2g^{k) 



(87) 



Owing to the angular momentum conservation, the pho- 
tons with Mz = do not scatter in the lowest order be- 
cause such photons cannot cause a direct transition from 
the lower to the upper state. They may cause such a 
transition provided it is accompanied by a simultaneous 
emission of a photon with = — 1 but this is a higher 
order process. Indeed, in the fourth order the scattering 
amplitude for Mz — 0, as seen from Eq. (|82bp . does not 
vanish. In comparison with the standard Breit-Wigner 
resonance formula, our r(w) is equal to half-width. As a 
result of the angular momentum conservation, the ampli- 
tude with AIz = 1 is resonant but the one with = —1 
is not. These two amplitudes correspond to two possible 
time orderings illustrated in Fig. [51 

The amplitude for the scattering of a photon off a two- 
level atom in the second order is obtained from Eq. ([7^ 
for the transition matrix 



where 



Am? — llP- — AmKiuj') — 4mif (tj) 



A(c.) 



r(c.) 



P 

7i'.g^(^) 

2\uj\ 



dkg^{k) 



(89) 

(90) 
(91) 



The standard resonance behavior can be seen after 
rewriting f{uj) in a different form. Disregarding the 
square of A{uj) + ir(w) and its product with since 
they are both of the fourth-order, we can approximately 



decompose the scattering amplitude (189)) into the follow- 
ing sum of simple fractions: 



2m — A(ti;) — cj — iT{u}) 

^ g^(^) 

2m — A(w) + uj — iT{Lj) 



(92) 



The first term is clearly resonant and the second is not. 
This expression agrees with the equal sign prescription — 
the width T{uj) enters with the same sign in both terms. 
Thus, the equal sign prescription, advocated in Ref. p^ . 
is appropriate for the photon scattering amplitude but 
not for the polarizability, as will be explained later. The 
scattering amplitude off a dipole atom has the same gen- 
eral form as for the two-level atom, so that our analysis 
applies also to this case. 



IX. LINEAR RESPONSE THEORY 

We shall use the quantum linear response theory (cf., 
for example, @) to relate the Feynman propagators to 
the important physical characteristics of two-level sys- 
tems: the spin susceptibility and the atomic polarizabil- 
ity. Linear response theory describes the reaction of a 
quantum system to a weak external perturbation. In the 
linear response theory, changes in an expectation values 
of observables are expressed in terms of retarded prop- 
agators. In our opinion most of the controversies in the 
treatment of damping resulted from the lack of a clear 
distinction between the scattering of photons (described 
by the S matrix and the Feynman propagators) and the 
time evolution of the expectation values (described by 
the solutions of the Heisenberg equations of motion and 
the retarded propagators). For the spin system, the ob- 
servables are the spin (or the magnetic moment) compo- 
nents '^'^cr'^. For the two-level atom, the observable is 
the atomic induced dipole represented (up to a constant 
factor) by the operator ^'V^j^'. The spin susceptibil- 
ity determines the response of the magnetic moment to 
the applied magnetic field and the polarizability deter- 
mines the response of the atom to the applied electro- 
magnetic field. These external perturbations are repre- 
sented by the terms '^'^ (T^-5^p{t) or '^^ (Tx'^5ip{t) added to 
the Hamiltonian. The change in the average spin value 
produced by such a perturbation is Q 

5{G\S,{t)\G) 

/OO 
dt'e{t - t'){G\[S^{t), Sj{t')]\G)Sipj{t'). (93) 
-OO 

The spin operators in this formula are in the Heisenberg 
picture S{t) = ^>'' {t)cr'^(t). Since the ground state |G) 
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is stationary with some energy Eq, we can write 

{G\m),S,{t')]\G) 

= (G|5,(0)e-'(^-^«)(*-*''5,(0)|G) 

- (G|S',(0)e*(^-^'^'(*-*''5,(0)|G). (94) 

Thus, the integral in Eq. ((93|) is a convolution and we 
can transform this equation to a simple algebraic relation 
between the Fourier transforms 



(95) 



The spin susceptibility Xiji'^) describes the response of 
the spin to a monochromatic external magnetic field and 
is defined by the Kubo formula [20] 



dti 



-iujt-\-et 



{G\[Sm,Sm\G), (96) 



where the damping factor e*^* guarantees that the applied 
field has been switched on adiabatically. 

The corresponding formula for the polarizability of a 
two-level atom reads 



5{G\S^(t)\G) 



dt'e{t - t'){G\[SAt), S^{t')]\G)5^{t'), (97) 



a{uj)^iA dte-^^'+''{G\[SM.S.,{t)]\G). (98) 



For a single two-level atom, the constant factor A is usu- 
ally (cf., for example, [ll,[lj|) given the value A = d'^/3H, 
where d measures the strength of the dipole transition. 
We have reversed the sign in the definition of a(w), 
as compared to the definition of the spin susceptibility 
([96|) . to be in agreement with the standard Kramers- 
Heisenberg-Dirac expression for the polarizability. 

The expectation value of the retarded commutator of 
the spin operators appearing in (j93p is directly related to 
the retarded photon propagator Gmj 

Gr^j {k,k',t-t')^ -ie{t - t'){G\ {t) , $j {t')]\G) . 

(99) 

Indeed, with the use of the Hcisenberg equations of mo- 
tion (j31bp for the field ^> and the canonical commutation 
relations (j32bp . we obtain 

= ~{d^ + k^)5{t - t')5,j5{k - k') 

- ig{k)g{k')e(t - t'){G\ [5. (i), S, {t')]\G) . (100) 

After the Fourier transformation with respect to t ~ t' , 
this relation becomes 

{kl~k^)(kl-k'^)gR,,{k,k',u;) 

= {kl - e)5,,5{k-k')+g{k)g{k')x^■J{Lo)■ (101) 

Analogous relations hold for the two-level atom. Thus, 
the susceptibility and the polarizability are simply pro- 
portional to the retarded photon propagator. 



A. Spectral representation 

One the advantages of using the methods of relativis- 
tic field theory is that the analytic properties of propa- 
gators became explicit. The retarded photon propaga- 
tor cannot he calculated by a direct application of the 
Feynman-Dyson perturbation theory. However, once we 
determine the photon Feynman propagator, the retarded 
propagator can be unambiguously reconstructed. In or- 
der to prove this assertion, we shall follow the same pro- 
cedure that in a relativistic quantum field theory leads 
to the Kallen-Lehmann representation. Starting from the 
general definition of the photon Feynman propagator we 
arrive at the following formula (cf., for example, Q): 

gF^J {k,k',t-t')^ ~l9{t - t')(G\ (fc, t)$, (fc', t')\G) 

- i9{t' - t)(G\^,{k',t')^j{k, t)\G) 



■--i6{t-t') / dMe-'^*^'-*'^A,j{M,k,k') 
Jo 

/>OC 

ie{t'-t) / dMe-'^^^*-'">Aj,{M,k\k), (102) 



where the spectral matrix Aij{M, k, k') is defined as fol- 
lows: 

A,j{M,k,k') = AUM, k',k) 



5{M -En + £;G)(G|$,(fc,0)|n)(n|$,(fc',0)|G) 

(103) 



{|n)} is any complete set of stationary states of the sys- 
tem and Eq is the energy of the ground state. Therefore, 
the Fourier transform QFij{k,k' ,kQ) of the propagator 
can be written in the form of a spectral representation 

Gpijik, k', ko) 

/A,(M,fe,.0_A-.(M,.,fcOX 
\ko- M + ie ko + M-ie J ^ ' 

Repeating the same procedure for the retarded propaga- 
tor defined in Eq. ((99|) . we obtain its spectral represen- 
tation 



gRij[k,k' ,ko) 



dM 



A,j{M,k,k') _ A*^{M,k,k') 
ko - M + ie ko + M + ie 



(105) 



The spectral matrices are the same and the only differ- 
ence is the change of the sign of the e-term in the second 
part. Since M > 0, this is equivalent to the replacement 
of ie in the denominators of the Feynman propagator by 
isgn(/co)e. This "epsilon rule" in the simplest case of a 
one-component propagator reduces to the standard rule 
of quantum linear response theory ^ : 

RegR{k,k',ko) ^RegFik,k',kQ), (106a) 
lmgR{k,k',ko) ^sgn.{ko)lmgF{k,k',ko). (106b) 
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These simple relations do not hold in general because the 
imaginary unit may appear not only together with e but 
also in other places. In particular, they do not hold for 
the spin system in the Cartesian basis. However, in the 
angular momentum basis the photon propagator is diag- 
onal. Therefore, we can treat each component separately 
and use the relations (|106p to obtain the components 
of susceptibility from the Feynman propagator (or more 
precisely from the scattering matrix) 



Rexa{cj) = ReTa{Lj), 



(107a) 
(107b) 



where the subscript a takes on the values ± and 0. Here, 
unlike in the formula ()85|1 for the scattering amplitude, in 
Eqs. (|107p the frequency cu takes on positive and negative 
values because the real function ip{t) describing an exter- 
nal perturbation must contain both positive and negative 
frequencies. This important difference must be kept in 
mind when discussing the relations between the scatter- 
ing and the linear response. 

The relations between the Feynman propagator and 
the retarded propagator in the simple form ()106p hold 
for the two-level atom and for the dipole atom. In the 
first cases there is just one scalar function from the very 
beginning. In the second case, due to the conservation of 
the three components of angular momentum, the propa- 
gator is proportional to dij so that it effectively reduces 
to just one function. 

We shall confirm the validity of the spectral represen- 
tations of the photon propagator by direct calculations 
in perturbation theory. 

Since our theory is invariant under time reversal, the 
retarded propagator will automatically satisfy the gen- 
eral crossing relation 



Gmjik, k', -ko) = Gnjiik, k', ko 



(108) 



This relation for a two-level atom implies that the atomic 
polarizability satisfied the condition a*{uj) — a{~uj). 
The significance of this condition has been emphasized 
in Ref. [7,]. 



B. Spin susceptibility 

The components of the photon transition matrix T{kQ) 
evaluated in Sec. IVIII gives the following formula for the 
spin susceptibility in the lowest order of perturbation the- 
ory: 



(109a) 



2m T w - 2 (A(w) ± i sgn(w)r(w)) ' 
xl^\co)^0. (109b) 



From the transition matrix (|82p we can obtain the spin 
susceptibility in the fourth order 



(2+4)/ N 



2(1-6) 



2mTuj -S -2{l-b) (A(w) ± i sgn(a;)r(w)) ' 

(UOa) 



(2+4). s „(2+4), ^ 



(110b) 



There is only one term for each transition but the op- 
posite sign prescription is still visible. The sign of the 
imaginary part in the denominator depends on the sign 

of LJ. 



C. Atomic polarizability 

In this Section we shall use the Hamiltonian (fTTjl for 
the two-level atom to calculate the photon propagator 
up to the fourth order of perturbation theory and use it 
to find the polarizability. All Feynman diagrams corre- 
sponding to the radiative corrections that will be taken 
into account in our calculation are shown in Fig. The 
frequency-dependent atomic polarizability q;(w) can be 
obtained from the functions T(ko) and T(fco) by chang- 
ing their imaginary parts according to the prescription 
dMl). 

The second-order self-energy part for the two-level 
atom is given in Eq. ((76|) . Substituting this expression 
into the formula (j65p and changing the sign of the imag- 
inary part, like in Eq. (|107p . we obtain the following ex- 
pression for the atomic polarizability 



ArnA 



Am? — up- — Am I A(cli) + i sgn(a;)f (w) 



(111) 



In order to see that this expression obeys the opposite 
sign prescription we could write it in the form of a spec- 
tral representation. However, following the treatment of 
the photon scattering amplitude, we shall convert this 
expression into simple fractions neglecting again higher- 
order terms 



A 



2m — A(w) — uj — i sgn(u;)r(a;) 
A 



2m — A(ti;) + lu — i sgn(w)r(u;) 



(112) 



Depending on the sign of lu, either the first or the second 
term is resonant. Therefore, if we only care about the 
important resonant terms, we can write this formula as 



a(2)(a;) 



A 



2m — A{uj) — LU — iT{Lu) 
A 



2m- A{lu) +LU + iT{Lu) 



(113) 
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In other words, this expression is a good approximation 
to the exact formula (jllip near both resonances, when 
bj K, ±2m. Thus, our expressions for the atomic polariz- 
abihty derived from the quantum hnear response theory 
agree with the opposite sign prescription, as advocated 
in Refs. [i,[ll[ii]. 

To extend this result to the fourth order of perturba- 
tion theory, we use the formula ()83|1 for the photon prop- 
agator. The resulting expression differs from the formula 
(imp for the polarizability in the second order only by 
the presence of the factors (1 — 



4m(l - b)A 



Ani^ — up- — 



4m(l - h) {k[uj) -I- i sgn(tj)f (w) 



(114) 



Our result is quite different from the formula derived in 
Ref. 0. It seems to us that this difference is due to the 
difhculties in systematically accounting in the standard 
treatment for all higher order corrections. In particular, 
Loudon and Barnett have not included all corrections to 
the ground state up to the fourth order and they have dis- 
regarded all level shifts. In our formulation, the method 
of Feynman diagrams guarantees an unambiguous deriva- 
tion of all corrections in any order of perturbation theory. 

Finally, we would like to emphasize that all those equal 
or opposite sign prescriptions, that are widely used in 
the semiphenomenological treatment, have some practi- 
cal limitations. They can be directly applied only to the 
spectral representations (I104p or (|105p . In general, as 
seen for example in Eq. (|114p . the expressions obtained 
directly in perturbation theory cannot be easily decom- 
posed into two parts because they are not given in the 
form of a spectral representation. Of course, we can al- 
ways find this representation, but the formulas are quite 
complicated and they hide the resonant character of the 
process. Still, having closed expressions we can always 
identify the analytic properties of aiu:) that correspond 
to these prescriptions. Namely, we can locate the posi- 
tions of the poles of a{io) in the complex uj plane to dis- 
cover that the pole near 2m lies in the upper half-plane 
while the the pole near —2m lies in lower half-plane. This 
property extends the opposite sign prescription to the 
general case. 



X. CONCLUSIONS 

We have shown that the methods of relativistic quan- 
tum field theory applied to a two-level (and also to a 
many-level) system interacting with the quantized elec- 
tromagnetic field lead to significant simplifications in the 
evaluation of various physical properties of the system. 
Owing to these simplifications we can easily go easily be- 
yond the lowest orders of perturbation theory. The dif- 
ference in complexity of the calculations performed with 



the use of the traditional approach and the new methods 
is enormous. For example, the interaction Hamiltonian 
for the spin system has six terms, so there arc 6* = 1296 
terms in the fourth order of the standard perturbation 
theory while in our approach we have only several Feyn- 
man diagrams to consider. It is true that for a partic- 
ular process many terms will not contribute but still a 
lot of terms must be taken into account. In addition to 
the simplifications in the calculations, we also gain new 
physical insights that stem from the connections that ex- 
ist in quantum field theory between different characteris- 
tics of the system. In particular, the connection between 
the photon scattering amplitude and the linear response 
functions of the two-level system to an applied electro- 
magnetic field is very useful. This connection is crucial 
to the understanding of the hotly debated relation be- 
tween the equal sign prescription and the opposite sign 
prescription in the description of damping. 
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APPENDIX A: MULTIPOLE EXPANSION 

1. Multipole expansion of the electromagnetic field 

The decomposition of the electric and magnetic field 
operators into the eigcnfunctions of the angular momen- 
tum can be written in the form 12 ll. 12211 



Eir) 



E 

JMX 



dk 



E 



(Al) 



B{r) 



B 



(A) 
IMk 



(r)4^i(fc)-fB7^;i,(r)cy;/(A:) 

(A2) 

where J > and A takes on two values (e, m) that dis- 
tinguish between the electric and magnetic multipoles. 
The operators Cj^^,(A:) and c^jMki^) annihilate and cre- 
ate photons with the energy fuv = Hck, the square of 
the total angular momentum ffJ{J + 1), the projection 
of the total angular momentum on the z axis HM , and 
the parity as determined by A. These operators obey the 
standard commutation relations 

CjM(fc)' 4'M'(fc')] = Sjj,SMM'Sxvd{k - k'). (A3) 

The energy operator of the electromagnetic field ex- 
pressed in terms of the creation and annihilation oper- 
ators is 



(A) 



H 



f 



JMX 



dktuvc^}^]ikyf>,ik). 



(A4) 
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The functions appearing in this decomposition (|A1[) and 
(|A2|) can be expressed in terms of the foUowing solutions 
TjMk{r) of the scalar Helmholtz equation: 



TjMk{r) 



7rJ(J + 1) 



]j{kr)YjM{n), (A5) 



where jj{kr) is the spherical Bessel function, YjM{n) is 
the spherical harmonic, and n is the unit vector in the r 



direction. The functions Bf'^j^{r) and E yi^^, (r ) can be 



expressed in terms of TjMk{r) as follows |2lLl22l|: 



(A6a) 
(A6b) 
(A6c) 
(A6d) 



B 
E 
B 



pie) 

■'JMk 

(e) 
JMk 
(m) 
JMk 
(m) 
JMk 



(r) = iVxLTjMkir), 

(r) = kLTjMkir). 

(r) = kLTjMkir), 

(r) = -iVxLTjMkir), 



where L — —ir x V is the angular momentum operator. 



2. Magnetic dipole coupling 

In order to extract the relevant part of B{r) that will 
contribute to the interaction Hamiltonian, we shall insert 
the expansion (|A2p into the formula (j3bp and obtain the 
following integrals (and their complex conjugate counter- 
parts): 

d^r p{r)B'flj^ir) = J d'r p{r)kLTjMkir), (A7a) 



d'rp{r)B'f^l{r) 



jd\p{r)VxLTjMk{r) 



(A7b) 



The integral in Eq. (jA7a|l vanishes because after the in- 
tegration by parts the angular momentum operator L 
acts on the spherically symmetric function p{r). How- 
ever, the integral in Eq. (|A7b|) might contribute. In 
order to explore this possibility, we use the identity 
iVX-L = — iLxV — 2V. Again, after the integration by 
parts we find that the term iixV does not contribute. 
We are left only with the gradient term Vp(r) = p'{r)n 
and we obtain 

j d'rpir)B^J^lir) = 2 J d\ p' {r)nTjMk{r). (A8) 

The integral over the angles will give a nonvanishing con- 
tribution only when the spherical harmonic YjMin) is a 
linear combination of the components of the vector n, i.e. 
for J = 1. Therefore, only the magnetic dipole compo- 
nent of the magnetic field contributes to the interaction 
Hamiltonian (j3bp . Using Eq. (|A8p we can rewrite the 
interaction Hamiltonian (I3bl) in the form 



Hi = -fiijj^axP- jd\^j^ dk ^j,{kr)p\r) 
X n \Yim {n)cM (k) + Y*^ {n)c\j {k)] . ( A9) 



To simplify the formulas we introduced the following no- 
tation: 



c+(fc) - c(™)(fc), co(fc) = ci7(fc), c_(fc) = cr\(fc)- 

(AlO) 

Using the explicit expressions for the spherical harmonics 

3 —X — iy 



(m) , 



>j=i,i(n) = 
y7=i,_i(n) = 
Yj=i,o{n) = 
and the formula 

j dO, riin. 



Stt r 
3 X — iy 



IT r 
3^z 
47r r ' 



47r , 



(Alia) 
(Allb) 
(Allc) 

(A12) 



we can perform the integration over the angles and write 
the interaction Hamiltonian (I5bl) in the form 



/•oo 

i// = V^V,t/>/ dkg{k) 
Jo 



c,(fc) + ct(fc) 



2k 



(A13) 



The annihilation and creation operators in the Cartesian 
basis arc built from the operators (jA10[) as follows 

c.(.)^^-M^, ct(.).^^I^l^, 

Cy[k) - , S(fe)-* ^ 

c,(fc) = co(A:), cl{k) = cl{k). (AM) 
They obey the standard commutation relations 

c,(fc),c](fc')] ^5,,5{k-k'). (A15) 

The formfactor g{k) is defined as 
pk 



9{k) 



d r p {r)ji (kr) 



(A16) 



The integral in this formula is proportional to the Fourier 
transform of p(r) 



d r p' {r)ji{kr) = —An / drp'{r)r ji(fcr) 



, , , . d f siaikr) cos(fcr) 
= A-K I drp{r)—\ — rrr^ " r ^ ^ 



dr \ fc2 



47r / drr p(r) smikr) 
Jo 

k I d^re''''''p{r) = kp{k). 



(A17) 
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3. Electric dipole coupling 

The calculation of the electric dipole Hamiltonian 
starts from the following formulas that are the counter- 
parts of Eqs. H]): 

+ \jd?r: {E\r) + B^{r)) : , (Al8a) 
H, = -e ( (fr'ip\r)ri}}{r)-E{r), (A18b) 



where is the atomic Hamiltonian. This time we will 
have four terms in the expansion of the electron field 
operators 



tp{r) 



(A19) 



where we used the Cartesian basis for the three wave 
functions of the degenerate upper energy level. The wave 
functions of the excited states and the ground state ap- 
pearing in this decomposition belong to the eigenvalues 
nie and uig of the atomic Hamiltonian. We continue to 
denote the energies by the letter m to stress the analogy 
with the relativistic QED. The value of the dipole mo- 
ment of the atomic transition between the ground state 
and an excited state, say the state described by xipe (r) / r, 
is 

d = e j tfr^ipe{r)xipg{r) = | J d^r ipeir)r(pg{r). 

(A20) 



Upon substituting the expansion (|A19P into 
Eqs. (|A18|1 . we obtain 

ffo = ip^'rhiP + i Jd\: {E\r) + B^r)) : , (A21a) 

Hi = -dip^'TXp- [ d^rK{r)n{n-E{r)) , (A21b) 



where n{r) = e(pe{r)r(pg{r)/d. Note, that according to 
Eq. (|A20p . K{r)/3 may be viewed as the normalized ra- 
dial distribution function of the dipole moment. It plays 
the same role as p{r) played in the description of the spin 
system. The diagonal matrix rh describes the energy lev- 
els and the matrices r describe the transitions between 
the ground state and the excited states 
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(A22) 



Substituting the expansion into multipoles (jAip of the 
electric field operator, we obtain two sets of integrals (and 
their complex conjugate counterparts) 



d\K(r)n{ri-E%,^(r)) 

J d^rK{r)n{in-VxLTjMkir)), 

d\n{r)n(n-E'-^l{r)) 

j d^r K,{r)n{n-LkTjMk{r)) . 



(A23a) 



(A23b) 



All integrals in the second set vanish because n-L = 0. 
The integrals in the first set can be simplified with the 
use of the relations in-'V xLTjMkir) = —L'^TjMk{r) = 
— J{J+l)TjMk{r)- We can again argue, as in our discus- 
sion of Eq. (jASp , that the only nonvanishing contribution 
in the formula (jA23a|l comes from the dipole component, 
when J — \. Thus, the interaction Hamiltonian becomes 



d^r^j dk]J —ji{kr)K{r) 



M 



X n 



YiM {n)dM (k) + Y*M {n)dlj (fc)! , (A24) 



where we have introduced again a simplified notation for 
the annihilation and creation operators 



d+{k) = do(fc) = c'llik), d^{k) = c'^'ik) 



(e), 



.(e) 



(A25) 



Note, that the electric dipole Hamiltonian (|A24p has the 
same general form as the magnetic dipole Hamiltonian 
(IA9p . Therefore, we may use the same methods to trans- 
form (|A24p to the form 



roo 

Hi = i^^Tiip / dk g{k) 
Jo 



d,{k)+d\{k) 



2k 



(A26) 



where 



gik) = -^^^ I d^rK{r)ji{kr) 



(A27) 
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APPENDIX B: RADIATIVE CORRECTIONS TO THE PHOTON PROPAGATOR 

1. Spin system 

The contributions to the photon sclf-encrgy in the fourth-order of perturbation theory corresponding to the diagrams 
(b), (c), (d), and (e) in Fig. [6] lead to the following integrals: 

/ ■7r'^T^WaSFiPa + ko)(^mSFiPO + ko + lo)(^bSFiPO + lo)(^nSFiPa)}DFmn{lo), (Bl) 
-oo ZTT 

T^'^ab\ko) = ( ( ■^^T:WaSFiPO + ko)abSFiPo)crmSFiPO + lo)'^nSF{po)}DFnmilo), (B2) 



-OO 

nOO 



'PlT\ko) ^ i Snis I ^TilaaSFipo + ko)abSF{po)crzSFiPo)} (B3) 

'PiT\ko) = -init [ ^^Ti'laaSFipo + ko)abSFipo)<7zSFiPo)} ■ (B4) 
J-oo 27r 

The contributions corresponding to the diagrams (f), (g), and (h) are the same as those corresponding to the diagrams 
(c), (d), and (e). Therefore, we shall calculate the contributions only from the diagrams (c), (d), and (e) and multiply 
them by 2. Again, as in the lowest order, the matrices Vab sue diagonal in the angular momentum basis. All Feynman 
integrals can be evaluated by the method of residues. We shall present the detailed calculation of the first integral. The 
three remaining contributions are even simpler to calculate and we shall give only the final results. The nonvanishing 
components of 7^^^''^(fco) in the angular momentum basis are V^''\kQ) and VQ'^''\k()) 

vf\ko)= rdkg^k) — L 







27r 27r - fc2 + ie 



^"^^'"i^- — TT — 7 ^1 — — n — rT~i n — ^+ — n — \ — r~ — n — r (^^'^) 

Po + kQ-(m-ie)az po + ka + Iq - [m - leja^ po + Iq - [m - lejcjz po-(m-te)az) 

dpo f°° dlo 1 1 1 1 1 



27r Iq — k? + ie po + kg — m + ie po + ko + Iq — m + ie po + Iq + m — ie po + m ~ ie^ 



pW(fco)= rdkg^k) — L 



-oo 27r 2tt l^-k^+ ie 

xTr<jcr+ ] ^ — cr^ ; — cr_ ; j —a^ ^ ^ — I (B5b) 

po + ko - [m~ ie)az po + kg + Iq - [m - lejcjz po + h - [m - ie)az po- (m- ie)az ) 

dpo [°° dh I 1 1 1 1 



27r Iq — k"^ + ie po + fco + m ~ ie pa + ka + Iq + m — ie po + Iq — m + ie po ~- m + ie' 
dpo f°° dlo 1 



27r - A;2 + ie 



1 1 11 

Tr <^ CTz ; — <7+- 



'po + ko - (m - ie)cr^ po + ko + Iq - {m - ie)(Jz po + Iq - {m - ie)(j^ pq - {m - ie)<Tz 

^ , 1 1 1 1 

Tr < ffz ; ; — ct- ; ; ; — f,- 



' Po + fco - ("^ - i()cz Po + ko + lo- (m- ie)az " Po + lo - {m - ie)az Pa - (m - ie)az 



(B5c) 



OO 



:-2 j dkg^ik) 



dpo [°° dh I 1 1 1 1 



27r 27r Zq — fc^ + ie pq + k^ — m + ie pq + ko + Iq + m — ie po + Iq + m — ie po — m + ie 

dpo dlo I 1 1 1 1 

^2tt Iq — k'^ + ie po + ko + m - ie po + ko + lo — m + ie po + lo — m + ie po + m — ie 



Note the advantage of using the angular momentum bases — in each case we are left with only one or two integrals. 
We dropped several terms that did not contribute because they contained the squares of the a± matrices. Finally, 
we replaced all az matrices by ±1 due to the presence of the projection operators Pg and Pg. The results of the 
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integrations are 



[Zm — Koj Jo k K + Zm — ko ~ le 

ri"'\ko)^rf\-ko), 



Iq k k + 2m {k + 2mY — k^ — ie' 
The contributions corresponding to the diagrams (c), (d), and (e) are 

2 f°°dk g^ik) 1 r°°dk / 1 



vi''\ko) 



2m — ko Jo k {k + 2m)^ (2m 

dk g^{k) 



1 rdk 



k + 2m — ko — ie fc + 2m 



1 



k k + 2m{k + 2mY - k^ 



vi''\ko) 



vi''\ko) 



2 5ms 



{2m-koY 
2mt 



, r[''Hko) = vi''H-ko), 7'r^(fco) = 



(B6a) 
(B6b) 
(B6c) 

(B7a) 
(B7b) 
(B7c) 
(B7d) 



(2m- fco)2^ 

Collecting all contributions in the second and fourth order, we obtain the following formulas for the three angular 
momentum components of the photon self-energy part: 

1 S 



vt+^\ko)^-2{l-h) 
Vr'\ko) = - 



2m =F ko {2m =F ko)^ 
dk g^ik) 1 



where 



b = 2 



k k + 2m (fc + 2m)2 — k^ — ie' 

dk g^{k) 1 f°° dkg^{k) 

T (fc + 2m)2 ' ~ 2(1-6) Jo P ■ 



(B8a) 
(B8b) 

(B9) 



2. Two- level atom 



The calculations of the fourth-order corrections to the photon propagator for the two-level atom are much simpler 
than those for the spin system. First, the propagator has only one component. Second, the contributions corresponding 
to the tadpole diagrams vanish. Third, there is no summation over three different states in the internal photon lines. 
The contribution corresponding to the diagram (b) in Fig. [6] is 



V^^''\ko)= / dkg^k) 



°° dpo [°° dlo 1 



X Tr (7, 



-oo 27r 2tt ll-k^+ ie 
1 1 



1 



Po + ko- (m - ie)cr^ po + ko + lo - [m - ie)(Jz po + lo - [m - ie)az po - [m - ie)a, 

dpo dlo I 1 1 1 1 



(BIO) 



f dkg^ r^jur^A 



Iq ~ k^ + ie po + ko + m — ie po + ko + lo — m + ie po + lo + rii — ie po — m + ie 



dkg^ik) 



°° dpo f°° dlo 1 



1 



1 



1 



1 



Q J _^ 2tt J _^ 2tt Iq — k"^ + ie po + ko — m + ie po + ko + lo + m - ie po + lo — m + ie po + m — ie 

(Bll) 



Similar integrals are obtained for the diagrams (c) and (e) in Fig. [SI The results of the integrations are 

2 f°^dk g^{k) 



v^^''Hko) = 

V^^-\ko) = 



Am^ — k"^ Jo k k + 2m ' 



1 



dk- 



8m2 



Am^-klJo (fc + 2m)2 [Am^~klYJo k k + 2m' 



'dk g^{k) 



r^^''\ko) = Sma 



16m 



4m2 - fc^ (4to2 - fc2)2 



(B12a) 
(B12b) 
(B12c) 
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The sum of the contributions from all diagrams in the second and fourth-order is 

73(2+4) (fco) ^ ^(2)(^^) ^ ^(4^)(fco) + 273(4c)(fco) + 2V^^^\k^) = -^y-^, (B13) 

where 

^ 1 /■°°dfc g2(fc)(fc + 4m) 

2m io /c (fc + 2to)2 ■ ^ 



3. Electric dipole atom 



For completeness, we also present the self-energy parts in the fourth order for the dipole atom. They are not much 
different from those for the two-level atom and we list here only the final results. There is one difference that is 
worth mentioning. Namely, the contributions corresponding to the diagrams (c) and (d) are now different from the 
contributions corresponding to the diagrams (f) and (g). This difference is due to the triple degeneracy of the excited 
energy level that breaks the symmetry between the |e) and \g) states that existed for the spin system and the two-level 
atom. Again, there are no contributions from tadpole diagrams. The self-energy parts in the fourth order are 

•pW^^ 2 f^dk g\k) 



2 


Ato2 - 


ftp 


2Am + 




Am2 — 


^-2 


2Am - 


ka 


Am2 — 


^-2 



K Ja '"kik + Amy 

where we used the values (|59p for the mass corrections. The total self-energy part in the second and fourth order is 
(disregarding the Kronecker symbols) 

p(2+^)(fco) = V^'\ko) + V^^'Hko) + V^'^Hko) + V^''\ko) + V^'f\k,) + V^^<^\k,) = (B16) 



where 



^ 1 f°"dk g^{k){k + 3Am) 

Am Jq k {k + Amy ' ^ ' 



[1] B. Schumacher, Phys. Rev. A 51, 2738 (1995). 

[2] E. T. Jaynes and F. W. Cummings, Proc. IEEE 51, 89 (1963). 

[3] R. Loudon, The Quantum Theory of Light (Oxford University Press, Oxford, 1973). 
[4] M. Gell-Mann and F. Low, Phys. Rev. 84, 350 (1951). 

[5] A. L. Fetter and J. D. Walecka, Quantum Theory of Many Particle Systems (McGraw-Hill, New York, 1971). 

[6] L. Allen and J. H. Eberly, Optical Resonance and Two-Level Atoms (Wiley-Interscience, New York, 1975), reprinted by 
(Dover Publications, New York, 1987). 

[7] R. Loudon and S. M. Barnett, J. Phys. B: At. Mol. Opt. Phys. 39, S555 (2006). 

[8] D. L. Andrews, S. Naguleswaran, and G. E. Stedman, Phys. Rev. A 57, 4925 (1998). 

[9] A. D. Buckingham and P. Fischer, Phys. Rev. A 61, 35801 (2000). 
[10] G. E. Stedman, S. Naguleswaran, D. L. Andrews, and L. C. Davila Romero, Phys. Rev. A 63, 47801 (2001). 
[11] A. D. Buckingham and P. Fischer, Phys. Rev. A 63, 47802 (2001). 

[12] D. L. Andrews, L. C. Davila Romero, and G. E. Stedman, Phys. Rev. A 67, 55801 (2003). 

[13] P. W. Milonni and R. W. Boyd, Phys. Rev. A 69, 23814 (2004). 

[14] P. R. Berman, R. W. Boyd and P. W. Milonni, Phys. Rev. A 74, 53816 (2006). 

[15] P. A. M. Dirac, Phys. Rev. 139 B, 684 (1965). 

[16] L Bialynicki-Birula and J. Piasecki, J. Math. Phys. 9, 1602 (1968). 

[17] F. J. Dyson, Phys. Rev. 75, 486 (1949). 

[18] J. D. Bjorken and S. D. Drell, Relativtstic Quantum Fields (McGraw-Hill, New York, 1965). 



24 



[19] C. Itzykson and J.-B. Zubcr, Quantum Field Theory (McGraw-Hill, New York, 1980). 

[20] R. Kubo, J. Phys. Soc. Japan 12, 570 (1957); R. Kubo, M. Toda and N. Hasliitsunic, Statistical Physics II, Nonequilibrium 

Statistical Mechanics (Springer, Berlin, 1985). 
[21] I. Bialynicki-Birula and Z. Bialynicka-Birula, Quantum Electrodynamics (Pergamon, Oxford, 1975), p. 172. 
[22] J. D. Jackson, Classical Electrodynamics (Wiley, New York, 1999, 3rd ed.), p. 431. 



